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Abstract 

The diffusion equation is extended by including spatiabtemporal memory in such a manner that 
the conservation of the concentration is maintained. The additional memory term gives rise to 
the formation of non-trivial stationary solutions. The steady state pattern in an infinite domain is 
driven by a competition between conventional particle current and a feedback current. We give a 



^ ■ general criteria for the existence of a non-trivial stationary state. The applicability of the model 

is tested in case of a strongly localized, time independent memory kernel. The resulting evolution 
C^ ■ 
H I equation is exactly solvable in arbitrary dimensions and the analytical solutions are compared with 

-h : numerical simulations. When the memory term offers an spatially decaying behavior, we find also 



the exact stationary solution in form of a screened potential. 
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I. INTRODUCTION 

There is an increasing effort in including feedback-couplings into the underlying evolution 
equations. Such memory-controlled effects should be a further unifying feature of complex 
physical |l|, |2| as well as biological systems 3| far from equilibrium. Whereas most of the 
papers are addressed to a purely time dependent but homogeneous memory, the present 
one gives an extension to spatially dependent processes. Such systems, however without 
memory effects, are widely discussed in physics j^ and biology [5]. Based on our recent 
studies of several evolution models ja, LD, |8|, |9| , where such non-Markovian memory effects 
had been enclosed, and on the analysis of chemical reactions [lO], we focus now on the 
diffusive behavior with memory couplings. These effects have an essential influence on the 
long time behavior, but they have also an impact on the dynamics in an intermediate time 
regime. As the consequence of the interplay between conventional diffusive transport and 
feedback we demonstrate the occurrence of non-trivial stationary solutions forming patterns. 
Our analysis should be grouped into the continuing interest in feedback couplings. A memory 
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behavior is well established in analyzing the freezing processes in under-cooled 
12| |. where the underlying mathematical representation is based on a projector 
formalism proposed by Mori [l3|. As the result of the projection procedure the irrelevant 
degrees of freedom, fluctuating on short time and short length scales, contribute to the 
time evolution of the relevant degrees of freedom as well as by instantaneous and by delay- 
controlled terms. Since the projection formalism is not restricted to selected systems, the 
modification of the evolution equation due to memory effects is rather general. Here the 
relevant quantity may be the concentration of certain species or the probability density for 
a particle. The crucial influence of memory effects can be illustrated by considering a single 
particle moving in a disordered environment. Due to the strong disorder each member of an 
equally prepared ensemble makes experience by its own surrounding, which is even modified 
by the random walk of the single particle itself. Such a behavior can be described by a non- 
linear Fokker-Planck equation with an additional memory term as had been demonstrated 
in [1J|. The analytical results, based on a one-loop renormalization group analysis, is in 
accordance with numerical simulations p and further analytical studies 7[. Whereas in 
that approach the memory effects are believed to be originated by the inherent non-linear 
interaction of the many particle system themselves, i. e. the time scale of the memory 



is determined by the relevant variable itself, there is a broad class of models which are 

.„...,„..,.,....Hya.......,„.„,.„,.„„.,„ 

|18| . Even in traffic models such memory effects has been discussed |19| . 



Recently j20|] memory effects in correlated anisotropic diffusion are studied in nanoporous 
crystalline solids. Likewise the effects of transport memory are discussed within the Fisher's 
equation 



21[, also applicable for bacterial population dynamics [22]. There appear a non- 
linear damping and traveling wave solutions 23|. The transport with memory, depending 
on the survivability of a population, is analyzed in 24 1. 

Regarding the large variety of systems with feedback couplings it seems to be worth to 
study simple models, which still include the dynamical features of evolution models as 
conservation of the relevant quantity p{f, t) and moreover, a time-delayed coupling. In the 
present paper we emphasize in an analytical solvable model with spatio-temporal feedback 
coupling. The generic behavior of the system may be changed by those additional delay 
effects. In particular, we discuss the transport behavior which is realized after a sufficient 
accumulation time and after cumulating particles within a spatial region. We show that the 
processes are affected by additional spatial and temporal correlations leading to non-trivial 
steady state solutions. 

II. MODEL 

In case the transport is realized after a spatial-temporal accumulation process the time 
evolution of the probability could also depend on the history of the sample to which it 
belongs, i.e. the changing rate of the probability should be influenced by the changing rate 
in the past. Thus the evolution of the probability p{f, t) has to be supplemented by memory 
terms. Such a term models, for example the way on which a seed probability at a certain time 
t' had been accumulated by a delayed transport mechanism, originated by the surrounded 
environment of the particle. In general, the changing rate of p at time t is also determined by 
the accumulation rate at a former time t' < t. In between, i.e. within the interval t = t — t', 
the particles are enriched while changing the probability at t'. Regardless that process the 
available amount of particles at time t is governed by an instantaneous transport term as 
well as by the changing rate at former times t'. Consequently the conventional diffusion 



equation should be modified according to J25| 

t oo 

dtp{r,t) = M{f,t]p,Vp) + dt' / d'^r'IC{f-r,t-t';p,Vp)C{f^,t';p,Vp). (1) 

~oo 

This equation is of convolution type. Whereas the operator A4 characterizes the instanta- 
neous and local processes, both operators /C and C are responsible for the delayed processes. 
In general the operators A1,/C and C may be also non-linear in p{r,t) and Vp(r', t). They 
have to be specified according to the physical situation in mind. In particular we show that 
the operators are restricted when p{f, t) is assumed to be conserved. 

A. Conservation 

In this section we specify the model, defined by Eq. (^, under the assumption that the 
basic quantity p(r, t) is conserved. This condition is fulfilled by 

dP{t) d ''+~ 



dt dt 



d\p{f, t) = . (2) 



To preserve p the instantaneous term Ai has to be related to a current, e.g. Ai oc V ■ j. 
Choosing natural boundary conditions, i.e. p = at the boundary, we find after Laplace 
transformation P{t) -^ P{z) 

zP{z)-Po = k{z)L{z) with k{z)= f d'^r}C{f,z), L{z) = f d\ C{r, z) , (3) 

with Pq = P{t = 0). For an arbitrary polynomial kernel K the conservation law is in general 
not fulfilled provided the operator C is simply defined hj C = —dtp{f, t) (the minus sign is 
only for convention). Making this ansatz we conclude from Eq. 1^ 

[zP{z)-Po][k{z) + l] = 

and consequently the conservation law is guaranteed. Physically, the assumption for C 
means, that we take into account a coupling of the rates, e.g. the evolution at the observation 
time t is directly coupled to the changing rate at t' < t. Processes in the past will be 
permanently "reevaluated" at present time. In doing so the memory kernel gives rise to a 
coupling of the time scales. In the vicinity of the upper limit t' ~ t the memory term reads 
/C[r, 0; p(r, 0)]dtp{f, t), i.e. a instantaneous change at the observation time t is coupled to the 
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value of p at the initial time t = 0. Therefore, the very past is related to the instantaneous 
value of p{r,t). In the opposite case, at the lower limit t' ~ 0, the change of the quantity 
p(r, t) near to the initial value df p{r, t' = 0) is directly coupled to the instantaneous value 
p{f, t) via the kernel. In such a manner the memory part represents a weighted coupling 
of the behavior at the initial time and the observation time. Due to that coupling of the 
rates the long-time behavior of the system will be modified. One reason for that could 
be that the moving particle is embedded into an environment of all the other particles of 
the system. Due to the mutual interaction, the particle is released after a certain waiting 
time t — t'. Especially in sufficiently complex diffusion- reaction systems the feedback and 
memory effects should be relevant. In such system additional degrees of freedom 
flexible macromolecules in melts or in concentrated solutions |26|1 , nematic elastomer 
in biology [18|]. In that context one is interested in the description of effiuent reprocessing 
plants in systems with closed water circulation. A special ecosystem of aerobic and anaerobic 
microorganisms is evolved in the clarifiers of such systems due to natural immigration or due 
to additional allowance. The living conditions of the microorganisms are mutually associated 
via the exchange of intermediate catabolic products. Each change of the concentration of 
one species will be stored in the food chain and effects the evolution of these species at 
a later time. Furthermore, the partial mixing in the clarifiers by convection and diffusion 
processes enlarges the effects over the whole system, so that the memory integral introduced 
in Eq. (^J) includes both, the time and the spatial coordinates. This special example may 
be extended also to other complex biological, chemical or engineering problems with various 
hidden degrees of freedom, which are able to influence the evolution of a selected component 
signiflcantly, for instance by biological interaction with other species via the food chain or 
via biological competition. Such effects which are partially observable, could contribute to 
the memory term. 

B. General properties 

Starting from Eq. (0) we can make some general statements for an arbitrary kernel 
/C(r, t;p, Vp) and the specialized term C{f,t;p,'Vp) = —dtp{r,t). As stressed before the 
last condition guarantees the conservation oi p{f,t). After Fourier transformation with re- 
spect to the spatial coordinate and Laplace transformation with respect to the time, we get 



from Eq. (^J 

p{k,z) = ^^ with b{k,z) = ^^, (4) 

z + D{k,z)k^ l + }C{k,z) 

where ]C{k, z) is the Fourier-Laplace transformed kerneL Eq. (j^ describes diffusion with a 
modified diffusion parameter D ^ D{f,t). Assuming a regular behavior of the kernel let us 
make the ansatz 

bik) 
]C{k,z) = ^^ + A{k,z) with lim A(A;, 2) = finite . (5) 

z 2^0 

Provided the kernel reveals a finite stationary value lim ]C{k,t) = b{k), then p{f,t) yields a 

t— >oo 

stationary solution, too. Inserting Eq. in Eq. (jlj) we obtain 



h{k) + Dk^ [h{k) + Dk^][z{l + K{k,z)) + h{k) + Dk^] 

(6) 

Summarizing the results we conclude that the model, following Eq. (Q) for the conserved 
quantity p{f,t), gives rise to a non-trivial stationary solution g{k), or equivalent g{r), if the 
Laplace transformed kernel /C(fc, z) satisfies the condition 

limz)C{k,z)^0. (7) 



z^O 



0. 



This result is a generalization of a previous one derived for a homogeneous system 

III. DIFFUSION WITH TIME INDEPENDENT MEMORY KERNEL 

In the following let us regard the evolution equation 

/oo pt 

d'^x' dt'K{f-f',t-t')dMf',t'). (8) 

-00 Jo 

As discussed in the previous section the conservation of the quantity p{f, t) is guaranteed 

for an arbitrary kernel K. Different to a conventional evolution equation, in the present 

one p(r, t) is related to p{f, t') with < t' < t. Moreover, Eq. (jS)) offers also a coupling 

between dtp{r,t) and dt>p{f,t'). The mixing of time scales leads, above all, to a substantial 

modification in the long time limit. 



A. Spatially localized memory 



As the simplest example we analyze a strictly spatially local, but time independent kernel 



Kif,t)=fx5ir), (9) 

where /i > is the strength of the memory. Note that in our approach spatial and temporal 
variables are decoupled. In case of a coupling between space and time variable one may 
observe different effects discussed elsewhere J28|. Inserting this kernel in Eq. (jHlwe get 



dtp{f, t) = DV'^p{f, t) — iJ.[p{f, t) — po{r)\ with po(0 = p(^5 1 = 0) 



(10) 



Such a memory term leads to a very strong coupling of the instantaneous value to the initial 
value because all previous time steps are weighted equally. There appears a memory induced 
feedback to the initial distribution as a driving force. Note that without the coupling to 
the initial value ^0(^)5 the particle would perform a diffusive motion and the probability 
density p(r, t) would decay on the time scale /i~^. Mathematically the inhomogeneous linear 
equation may be solved analytically for arbitrary initial conditions and in any dimension: 



p(f, t) = e''^* / G(f - f', t)po(r ') cf^r' + ii e"^* 



00 rt 
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G{f- f', t - t') e^*' poif')dt' d'^r' . 

(11) 
From this solution one can verify that a non- negative initial distribution po(^) ^ leads 

always to p{f,t) > provided // > 0. Eq. (fTTj) enables us to calculate the second moment 

s{t). It results 



s{t) = / r^p{f,t)d'^r 



2dD (1-e-^*) 



Po{f) d^r + r ^ Po{r) d'^r . 



(12) 



Notice that for /i -^ the last equation shows the conventional diffusive behavior. 

Let us illustrate the behavior in the one dimensional case with the initial condition po{x) 

Po6{x). The exact solution reads 



p{x,t) 



V4:TrDt 4 



erf 



(- 



X 



+ ^ i e- 



erf 



X 



V4Dt 



VViDt 
+ y/JIt] -sgn(-x) > , 



/it - sgn(x) 



(13) 



where erf (x) is the error function and n = ^ [ijD. The first part of Eq. ()13|) is the solution 
of the homogeneous equation showing a temporal decay with the time constant /x~^. For 



different strength of the feedback couphng /i = 0.5, 1, 3 the solution p{x,t) is presented 
in Fig. (P). The stronger the feedback strength n the more pronounced are the deviations 
from the pure diffusive behavior. In particular, the system develops a non-trivial stationary 
solution whenever /i > 0. Such a stationary solution is due the permanent coupling to the 
initial distribution. We find 

lim p{x, t) = p,{x) = ^e-'^l"! . (14) 

where Ps{x) is also included into Fig. (^. Notice that in the limiting case k —>■ oo the 

stationary solution tends to the initial distribution. Obviously our feedback offers a kind 

of shape memory. The stronger the memory is the stronger is the recognition of the initial 

state. 

The result can be generalized to arbitrary dimensions d. For the stationary solution Psir) 

we obtain with r =| f^ | 

PO d+2 d-2 . ^ 

P^y^) = f27rY/^ '^ ' ^ ' K^{Kr) . (15) 

The modified Bessel function K^{x) can be expressed by standard functions for d = 1,3 J29|. 
In the two dimensional case the modified Bessel function offers a logarithmic behavior. The 
evolution equation for p{f, t) is integrable with 

"^■'^^ (AnDtY/^^ +^^° I [A7TD{t-f)Y/^ ■ ^ ^ 

The last integral can be calculated exactly for d = 1 with the result given in Eq. ()13p. 
Asymptotically we obtain 



i.2 



The analytical results are confirmed by numerical simulations. 



B. Gaussian and exponential initial conditions 

For an arbitrary initial condition pQ{r) the stationary solution can be directly calculated 
by solving equation V^Ps(r) = K^[ps(r) — pQ{r) ] with k^ = ^/ D. It results in 



Ps{^ = -^^^^ / rfV , J''^ ;^ K._^{n I f - r- I) . (18) 



In this subsection we study different realizations of the initial conditions. In case of a 
Gaussian initial distribution po{r) = poe~'^^ the results for different strength of memory 
/i for both p{x,t) as well as Ps{x) are shown in Fig. (j2)). When the initial condition is an 
exponentially decaying function po{f) = pQe~^^^, we find solutions presented in Fig. Q. 
Summarizing this section we have shown, that a simple feedback coupling lead to stable 
patterns which do not appear under a conventional diffusion term. 

IV. FADING DIFFUSION 

Now we regard as a further realization of Eq. (0) with a conserved p{f, t) a spatially 
decaying memory term by choosing the kernel in the form 

/C(f,t) =/iexp(-A I f I). (19) 

Performing Laplace transformation we get 

Inserting that result into Eq. © it results a non-trivial stationary solution, compare also 



Eq. 



Ps{k) = — pom!^>^_^. (20) 



In the hydrodynamic limit k —>■ one finds the Fourier transformed stationary solution 



P.(n = |5^/<'^-'^^fLA-^(M?--|) with .= y^. (21) 



Using the asymptotic expansion for the modified Bessel functions J29|, one concludes that 
the stationary solution offers a screened Coulomb potential decaying on a length scale k~^. 
Note that the screening length is determined by the memory strength /i as well as by the 
characteristic length A~^ of the spatial memory. The screening length increases when the 
memory strength /i decreases. Let us stress that the non-trivial steady state is originated 
by the competition between diffusion and memory and it is not driven by spatial boundary 
conditions. In that case there appears an additional length scale, the size of the sample L, 
which gives rise to further effects not discussed now. 



V. CONCLUSIONS 

In this paper we have extended the conventional diffusive transport by including non- 
Markovian memory terms within the frame of evolution equations. The additional terms 
are chosen in such a manner that the relevant variable p{f, t) is normalized. By this require- 
ment the form of the memory term is restricted to a class where the changing rate of p at 
the observation time t is coupled to the changing rate of p at a former time t'. Insofar the 
delay effect offers a long time memory due to < t' < t. Further the model exhibits also 
a long range memory because the kernel /C depends on f—f'. Here we have considered 
several realizations of the memory kernel with different behavior. As a new result there 
occur non-trivial stationary solutions forming patterns. Such patterns are absent for con- 
ventional diffusion in an infinite domain without boundary conditions. The reason for such 
a new behavior is by means of the above mentioned coupling of the rates. This time accu- 
mulation is further accompanied by an additional spatial accumulation, the effect of which 
is comparable to the effect a long-range interaction forces and consequently the results are 
basically independent of the spatial dimensions in according to scaling arguments. These 
many-body effects are shown to change the asymptotic behavior drastically. Due to the 
feedback-coupling of a particle to its environment, a subsequent particle, undergoing a dif- 
fusive motion, gains information from a modified environment. As the result there appears 
patterns controlled by the memory strength. 
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(a) Stationary solution 



(b) ^ = 0.5 




(c) M = 1 



(d) M - 3 



FIG. 1: Stationary solution Ps{x) and p{x, t) in d = 1 with D = 1 and pq = 1 for different strength 
of memory fi. 
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(a) Stationary solution 




(b) ^i = 0.5 





(c) M = 1 



(d) /. = 5 



FIG. 2: Stationary solution Ps{x) and p{x,t) for the Gaussian initial distribution with D = 1, d 
1, Po = 1, X = 1 and different jj,. 
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(a) Stationary solution 



(b) ^i = 0.5 




(c) M - 1 



(d) M = 5 



FIG. 3: Stationary solution Ps{x) and p{x,t) for the exponential initial distribution with D 
1, d = 1, po = I, X = 1 and different /x. 
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